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ABSTRACT 

We present optical and near-infrared images and spectra of ESO248-G10, the host of 
the giant radio galaxy MRC B0319-454. From near-infrared colours, the active nucleus 
is shown to be reddened by hot dust emission or dust extinction. Star forming regions 
are identified beyond a radius of Sarcsecs (8.1 kpc) where hot gas is prevalent. The 
optical spectrum shows gas rotation at speeds of up to ~ 350 km s"^ out to > 32 kpc 
along the radio axis. A model is proposed in which the giant elliptical is triaxial with 
the radio axis along the short axis, and the figure rotation is around the long axis. From 
the model, the angles of the principal axis are if) = 34°, 6 = 65°, and = 19° with axis 
ratios of B/A = 0.75 and C/A = 0.69. A second velocity feature from 5 to 15arcsecs to 
the north-east of the nucleus is proposed to be a merging gas-rich galaxy inducing star 
formation while settling into an orbit about the figure rotation axis. 

Subject headings: galaxies:active — galaxies:elliptical and lenticular, cD — galax- 
ies:individual:ESO248G10 — galaxies:kinematics and dynamics — galaxies:nuclei — 
infrared : galaxies 
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1. INTRODUCTION 

Radio galaxies are frequently found to have giant 
ellipticals as their hosts (Carballo et al. 1998 and 
references therein). In elliptical hosts with a dust 
lane, the radio jets are typically aligned perpendicular 
to the dust lane (Kotanyi & Ekers 1979; MoUenhoff, 
Hummel & Bender 1992). The sense of rotation of 
the dust lane about the radio jet axis can be parallel, 
anti-parallel or perpendicular to the stellar rotation 
in a stable system. This implies an external origin 
for the dust. The common occurrence of dust lanes 
in elliptical galaxies suggests that infalling dust has 
been trapped in stable closed orbits. It is not surpris- 
ing, therefore, that Heckman et al. (1986) found that 
50—75% of radio galaxies show evidence of ongoing 
or past mergers/interactions. The orientation of the 
radio axis and rotation of a dust lane compared to the 
stellar rotation are indicative of the merger history of 
the galaxy. 

Many elliptical galaxies show evidence of complex 
kinematics. Warped dust lanes, radio jets misaligned 
with the major and minor isophotal axes, and tum- 
bling around an axis perpendicular to the radio axis 
are all indicative of triaxiality in elliptical galaxies. 
While a warped dust lane could indicate that the 
dust has not yet settled into a preferred plane, van 
Albada, Kotanyi & Schwarzschild (1982) have mod- 
elled warped dust lanes as stable orbits in tumbling 
triaxial systems. Binney (1978) first proposed that 
elliptical galaxies might be triaxial. Since then it has 
come to be accepted that most, if not all, ellipticals 
are triaxial and many also show evidence of tumbling. 
Unfortunately, the intrinsic shape and orientation of 
such a system cannot be determined from surface pho- 
tometry alone. While the apparent axis ratios can 
put limits on the lengths of the intrinsic axes, stel- 
lar and ionised gas kinematics together with the dust 
lane orientation are essential for the determination of 
a triaxial shape, as they trace the gravitational po- 
tential. 

MRC B0319— 454 is a giant radio galaxy with an 
FR II (Fanaroff & Riley 1974) morphology. It was 
identified by Jones (1989) with the Molonglo Observa- 
tory Synthesis Telescope (MOST). The double-lobed 
radio structure has a projected separation of 2.51 Mpc 
(Ho = 50kms~-^ Mpc~^), with quasi-continuous jets 
at a position angle of 51° and radio knots mapped 
by Saripalli, Subrahmanyan & Hunstead (1994). The 
host galaxy, ESO248-G10, has been identified with a 



giant elliptical galaxy at a redshift of 0.0625, giving a 
distance of 337 Mpc. At this distance 1 arcsec corre- 
sponds to 1.63 kpc. With a major axis of ^ 100 kpc 
at the R = 24.5 mag arcsec"^ contour, it is an ideal 
candidate to investigate the host galaxy properties of 
a giant radio galaxy. 

The optical images of Saripalli et al. (1994) show 
ESO248-G10 to have a warped dust lane. The radio 
axis is misaligned with the optical isophotal minor 
and major axes by 34° and 56° respectively, hinting 
that the galaxy may be triaxial. In this paper we 
attempt to deduce the true orientation of the princi- 
pal axes to the apparent isophotal axes from the dust 
lane morphology and ionised gas dynamics. Further- 
more, colour analysis is used to find direct evidence 
of a merger or interaction. Section || describes the 
observations and reduction, with the results in sec- 
tion ^ A discussion of the dynamics and orientation 
of ESO248-G10 is in section |. 

2. OBSERVATIONS AND REDUCTION 
2.1. Images 

The near-infrared bands of J, H, and Kn are ideal 
for studying the circumnuclear regions of ES0248- 
GIO as infrared wavelengths can penetrate the en- 
shrouding dust. Mosaiced images centred on the nu- 
cleus of ESO248-G10 were taken in and 
i^n-bands (see Table |l|). The Kn image was obtained 
with IRIS (Allen 1992) at the f/15 cassegrain focus of 
the 3.9 m Anglo- Australian Telescope (AAT) at Sid- 
ing Spring Observatory on 1994 October 22. Pixels 
were 0.6 arcsec on the sky, with 128 x 128 pixels on the 
HgCdTe array. The mosaic was comprised of three, 
five-point patterns centred on the Kn point source 
with offsets of between 10 and 15 arcsec for the sur- 
rounding frames. The standard star HD20223 was 
observed immediately after the image frames. Sky 
frames were offset by between 80 and 190 arcsec from 
the galaxy centre. 

J and H images were taken with CASPIR (McGre- 
gor et al. 1994) on the MSSSO 2.3 m telescope at 
Siding Spring Observatory on 1998 October 1-2 and 
1999 February 4 respectively. The 256 x 256 InSb ar- 
ray was used in FAST mode with 0.5 arcsec pixels. A 
nine-point pattern of images made up both the mo- 
saics, with offsets of 8 arcsec for each frame surround- 
ing the nucleus. Sky frames were offset by 300 arcsec 
from the nucleus. The standard stars HD20223 and 
HD15911 were observed immediately before or after 
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Table 1: Log of Observations. 









Av. Seeing 




Integration Time 


UT Date 


Type 


A 


(arcsec) 


Photometric 


(b) 


1994 Oct 22 


image 


Kn 


1.4 


Yes 


750 


1999 Feb 4 


image 


H 


1.3 


Yes 


1000 


1998 Oct 1-2 


image 


J 


1.5 


Yes 


1740 


1995 Feb 7 


spectrum 


Opt. 


1.25 


Yes 


2000 


1995 Oct 31 


spectrum 


HK 


2.5 


No 


11640 



the image frames. 

Data reduction and image processing was done 
with four software packages: iraf (Tody 1993), FI- 
GARO (Shortridge 1993), karma (Gooch 1996), and 
MIRIAD (Sault, Teuben & Wright 1995). The /Cn-band 
IRIS image was reduced with the IRIS-specific tasks 
in FIGARO, while the J and if images used the CASPIR 
reduction package in iraf. In each case the frames 
were flat fielded using dome flats and then sky- 
subtracted. All three wavebands were cleaned of bad 
pixels with FIGARO clean and mosaiced with IRIS- 
MOS. At least four stars in the field were used to regis- 
ter the frames when mosaicing the J and Kn images. 
The H image, however, was taken with tip-tilt cor- 
rection and blind mosaicing was sufficient. All three 
images were taken in photometric conditions and were 
magnitude calibrated using the standard star obser- 
vations. 

2.1.1. Alignment of Images 

The J, H, and Kn images were aligned relative to 
each other by the following method. The Kn image 
was resampled from 0.6 arcsec pixel"^ to 0.5 arcsec pixel 
with IRAF MAGNIFY to match the J and 77 images. By 
cross-referencing six stars on the digitized sky survey 
(DSSI, blue) image in karma koord, coordinates 
were attached to the J mosaic accurate to 0.1 pixel 
rms. Coordinates were transfered to the H and Kn 
images in the same way but the Kn image was smaller 
so that fewer alignment stars gave a larger scatter. It 
was then necessary to shift the H and Kn images with 
IRAF IMSHIFT to align the relative pixels to the J im- 
age to < 0.2 pixel. J — H and H — Kn magnitude 
images were then made in MIRIAD by subsetting an 
identical region from each of the images, then con- 
verting to magnitude images and subtracting. A co- 
ordinate system was applied to H — Kn and J — H 
images by copying the fits header parameters from 
the parent images. 



The B and R CCD images were those used by Sari- 
palh et al. (1994). They were taken at the AAT, 
and are measured to have a seeing disk FWHM of 
2.65 arcsec and 1.8 arcsec for B and R respectively. 
Coordinates were transferred to both images by com- 
parison with the DSS image. The resulting images 
align with the J image to < 0.2 arcsec. The R image 
was convolved to have to same resolution as the B 
image, before a B — R image was made with the same 
relative astrometric accuracy. 

2.2. Spectra 

2.2.1. Optical 

Optical spectra were taken for us as service ob- 
servations at the AAT on 1995 February 7 (see Ta- 
ble 0). The RGO spectrograph was used with the 
25cm camera, 270R grating, and Tektronix Ik x Ik 
CCD. The slit was 192.5 arcsec long with a spatial 
scale of 0.76 arcsec pixel" ^; it was aligned at a posi- 
tion angle of 40°, rather than the requested 51° (ra- 
dio jet axis p. a.) and consequently is not along any 
physically significant direction. Sky subtraction and 
cosmic ray removal was done in FIGARO. Wavelength 
calibration was done with CuAr lamp spectra and the 
measured spectral resolution is 7.5A FWHM. The re- 
sultant spectral image shows significant rotation in 
the stronger emission lines [Nii]A6584 and Ha, which 
can be seen in Figure 0(a)&(b). Spectra were there- 
fore extracted along the slit to measure the velocities 
either side of the nucleus. To improve the signal-to- 
noise ratio away from the nucleus an increasingly wide 
extraction window was used. A 3-pixel spectrum was 
extracted at the centre, then spectra with 3, 4, 4, 5, 
5, and 5 pixels were extracted consecutively to the 
north-east, and 3, 4, 4, and 5 pixels similarly to the 
south-west. Two more spectra of 9 pixels each were 
centred 19 pixels north-east and 18 pixels south-west 
of the centre. Each of the resulting spectra was di- 
vided by that of the standard star HD26169. From 
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each of these spectra, the [Nli]A6584 and Ha hnes 
were deblended with IRAF ngaussfit by constraining 
the relative flux, position and width of the [Nii]A6584 
& 6548 hnes. Other hnes were measured with iRAF 
SPLOT. Figure |l|(c) shows the rotation curve formed 
from the [Nii]A6584 and Ha hne velocities measured 
from each of the extracted spectra. Each point on the 
curve is plotted at the center of the extracted window 
that formed each of the 13 spectra. It is notewor- 
thy that the lines could be measured to a radius of 
~ 20arcsec to the north-east, but only ~ 14 arcsec 
to the south-west. The S/N in the Mgib line was not 
sufficient to measure a stellar rotation curve. 

Optical extinction is typically measured from the 
Ha/H/3 ratio. The H/3 line is, however, not measur- 
able with sufRcient accuracy to give a meaningful es- 
timate of extinction. Furthermore, as the nucleus is 
below t he h ighest obscuration of the dust lane (see 
section extinction correction was not applied 

to the emission line measurements. The errors in 
the [Nii]A6584/Ha line ratios due to extinction are 
negligible as the wavelengths are close. Figure |l](d) 
shows the [Nii]A6584/Ha ratios for flve of the spectra. 
These ratios will be discussed further in section 4.2.2| . 



2.2.2. Near- infrared 

An HK echelle spectrum was obtained with IRIS 
on the AAT on 1995 October 31. The wavelength 
coverage is from 1.45 to 2.5 /im, with a resolution of 
400. The slit was 13 arcsec long, with a spatial scale 
of 0.79 arcsec pixel" ^, and oriented east- west. The 
galaxy nucleus was positioned 3.5 arcsec either side 
of the centre of the slit in alternate frames. Typical 
IRIS reduction involves subtraction of pairs of frames 
to remove sky and background galaxy. While this 
method was used here, it was discovered later (see 
section 4.2.2 ) that beyond 4 arcsec from the nucleus 
along the optical slit angle, the galaxy has starburst 
processes not found at the nucleus. Starbursts re- 
gions may also be present in the direction of the HK 
slit, at a radius of 7 arcsec from the nucleus, in the 
region used for sky subtraction. There is therefore 
a chance that the sky background spectra that were 
subtracted may have starburst lines not found at the 
nucleus. However, contribution from these possible 
lines is expected to be small, as the nucleus is so much 
brighter than the background galaxy. This possible 
contamination is taken into account in evaluation of 
the spectrum. 

The echelle frames were flat-fielded and cleaned in 



FIGARO. After coadding, spectra were extracted from 
the central 4 pixels. Two standard stars bsl006 and 
HD21473 were used to flux calibrate the spectrum 
with IRISFLUX. As the observing conditions were not 
all photometric, this flux calibration is only relative. 
Wavelength calibration was done with argon lamp 
spectra and wavelengths were corrected to a helio- 
centric frame. 

3. RESULTS 
3.1. Images 

In Figure ||, the R image shows a dust lane cross- 
ing the front of the galaxy to the north-east of the 
nucleus. This suggests the galaxy is tilted away from 
the sky plane to the north-east. The nucleus is there- 
fore positioned below the strongest band of extinc- 
tion. The near infra-red images are much more cen- 
trally condensed, with the .R'n— band image having an 
unresolved core coincident, within errors, with the ra- 
dio core reported by Saripalli et al. (1994). The Jim- 
age shows a notch-shaped indentation approximately 
4 arcsec to the east of the nucleus, which is not ap- 
parent in the H and Kn images. Overlaid J — H and 
H — Kn extinction images formed from the division 
of the J, H , and Kn images are shown in Figure ||(d). 
As expected, the strongest reddening in the H — Kn 
image is above the nucleus, along the dust lane. 

Aperture photometry results are given in Table I 
Apertures were chosen to extend well past the north- 
east edge of the dust lane, with measurements made 
using the IRAF RADPROF program. All apertures 
are centred on the ifn— band peak, which coincides, 
within errors, with the peaks in J— and i/— bands. 
As a background aperture is required for photome- 
try, and the three images are different sizes, an an- 
nulus was chosen that fits within each image and 
avoids any stars. With inner and outer radii of 36 
and 41 arcsec, this background aperture avoids the 
main galaxy halo, but a small contribution from the 
galaxy cannot be ruled out. At large apertures, this 
introduces an estimated systematic error of up to 
0.1 magnitudes. In smaller apertures, however, the 
error is dominated by standard star photometry er- 
rors, which are 0.02 magnitudes. 

Aperture photometry in i^Tn-band gave the ab- 
solute magnitude in a 1.5 arcsec radius aperture as 
M/f = —25.9 ± 0.5. The canonical luminosity criteria 
for optical QSOs is Mk — —25.8 (Surace & Sanders 
1999), while typical Seyferts only have Mb = -20.2. 
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MRC B03 19-454 [Nil] >.6584 & Ha Rotation Curve 
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Fig. 1.— (a)Optical Spectral Image including the [Nii]A6584 and Ha lines. The inner pixels clearly show fast 
rotation of ~ lOOkms"^ arcsec^^. (b) Same as (a) with greyscales chosen to highlight the extent of the rotating 
gas out to almost 20arsecs. (c) Rotation Curves for [Nii] (crosses) and Ha (stars) along p. a. of 40°. See text for 
explanation of the points, (d) Change in the relative strength of [Nll]A6584 and Ha with position. Each spectrum 
is for the region marked by the corresponding line on th^ rotation curve and has a width of ~ 182 A. 




Fig. 2. — (a)i?, {h)Kn, and (c) J band images, (d) H — Kn extinction image, with overlaid contours of the J — H 
extinction image; darker regions indicate greater reddening. Contour levels are 0.86-0.98 in 0.02 steps. 
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Table 2: Aperture photometry centred on the IR point 
source. Errors in the colours range from 0.02 — 0.1 mag 
from smallest to largest apertures. 















Radius 


J 


H 


Kn 


J-H 


H - Kn 


(arcsec) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


15.24 


14.30 


12.66 


0.94 


1.65 


2 


14.24 


13.35 


11.47 


0.89 


1.89 


3 


13.82 


12.96 


10.94 


0.86 


2.02 


4 


13.54 


12.71 


10.64 


0.83 


2.07 


5 


13.35 


12.54 


10.45 


0.81 


2.09 


6 


13.19 


12.40 


10.30 


0.79 


2.10 


7 


13.07 


12.30 


10.18 


0.77 


2.11 


8 


12.96 


12.21 


10.09 


0.76 


2.12 


9 


12.88 


12.13 


10.01 


0.74 


2.13 


10 


12.80 


12.07 


9.94 


0.73 


2.13 



Therefore, the infrared hght from ESO248-G10 is suf- 
ficient to class it with the QSOs. 

3.2. Spectra and Rotation Curves 

3.2.1. Emission Lines and Ratios 

Figure || shows the optical spectrum at the nu- 
cleus. The emission lines of [N ii]AA6548, 6584, Ha, 
[O iii]AA4959, 5007, [O i] A6300, and [S ii] AA6716, 6731, 
and the absorption lines of NalD and Mgib are all 
marked. 

In Figure U, the near-infrared HK spectrum has 
some absorption lines, but no emission lines are ev- 
ident. Significant lines that lie within the spectral 
range at this redshift include two hydrogen recom- 
bination lines (Paa and Br7), several [Fell], and H2 
lines. Hydrogen recombination lines are expected to 
be stronger from starbursts than AGNs due to the 
Hii regions surrounding young hot OB stars (Hill et 
al. 1999). This may suggest that there is not a strong 
contribution from starburst regions within a radius of 
1.6 arcsec (the offset of the extracted spectrum) from 
the nucleus. The lack of [Fell] and H2 lines is less 
indicative as their origins can be starburst or AGN. 
For starburst regions, [Fe 11] lines can be excited in the 
cooling tails of supernova remnants, or directly pho- 
toionized by young OB stars. Alternatively, X-rays 
from the nuclear source may photoionize the NLR 
clouds, or radio jets may shock-excite iron. In sec- 



tion 4.2.2 , we show that the optical line ratios in the 
nuclear regions provide strong evidence for AGN ac- 
tivity while further out, starburst processes are found. 



The lack of [Fe 11] lines in the nuclear spectrum indi- 
cates that the AGN processes are not exciting iron 
in this galaxy. There may be a lack of iron within 
1.6 arcsec of the nucleus, or alternatively, this could 
point to the starburst processes as more likely to be 
the excitation mechanism for iron in active galaxies. 

3.2.2. Rotation Curves 

The ionised gas rotation curve in Figure |l|(c) clearly 
shows fast rotation along the slit direction, which has 
a position angle 11° from the radio axis. The ac- 
tual rotation axis is unknown. The south-west side 
is approaching and the north-east (dust lane) side is 
receding. As the distribution of emission lines is not 
symmetric about the nucleus (see section [2.2.1 ), the 
ionised gas extends further to the north-east than the 
south west. The rotation curve shows a sudden de- 
crease in velocity from 5 ± 1 arcsec to ~ 15 arcsec to 
the north-east. A comparison of the features in Fig- 
ure H shows that the dust lane extends to a maximum 
offset of ~ 4 arcsec from the nucleus, which is before 
the decrease in velocity and is therefore unlikely to 
be related. Similarly, Nicholson, Bland-Hawthorn & 
Taylor (1992) found that the dust lane obscuration 
has only a minor effect on an Ha image of the dust 
lane elliptical galaxy, Centaurus A. The slowest veloc- 
ity point in the 'dip' does not have a corresponding 
feature in the ESO248-G10 images. 

It is noteworthy that the rotation curve shows no 
sign of a turnover out to 20 arcsec or 3.3 effective 
radii. Extended rotation curves have been measured 
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Fig. 3. — The optical spectrum of ESO248-G10 at the nucleus. Relative values only are shown on the intensity 
scale as the spectrum is not flux calibrated. 
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1.6 1.62 1.64 1.66 1.68 

Wavelength microns 




Wavelength microns 



Fig. 4. — Two regions of the near-infrared spectrum are shown, shifted to the rest frame. The expected positions 
of [Fe ii], ^^CO, and H2 hnes are marked. While some absorption hnes are evident, none of the expected emission 
lines are clearly detected. Typically, the strongest lines would be the [Fe 11] lines at 1.644/im and 1.664/im and the 
H2I — S(l) line. The wavelength calibration is accurate to 0.003 /xm. 
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Fig. 5. Graph of the Ha and [Nii]A6584 average rotation curve and the profiles across the R— and J5— band 
images along the slit. This is then overlaid on the iJ— image such that the graph radius-axis aligns with the slit 
position angle through the nucleus. 
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for relatively few ellipticals (e.g. NGC 3115, Capac- 
cioli et al. 1993) and in each case, have been found 
to remain flat out to 3 — 4 effective radii. 

4. DISCUSSION 

4.1. Galaxy Rotation and Orientation 

4.. 1.1. Orientation of the Principal Axes 

The emission line gas shows fast rotation along an 
apparent axis close to the radio axis (Figure |^(c)). 
Elliptical galaxies with rings or disks rotating about 
an axis perpendicular to the symmetry plane are not 
rare. A survey of 40 early- type galaxies with dust 
lanes by Hawarden et al. (1981), found that po- 
lar disk/ring galaxies comprised the largest subgroup, 
suggesting that they are likely to be stable systems. 
Such a polar gas disk is only stable in a triaxial sys- 
tem (Steiman-Cameron & Durisen 1982; van Albada 
et al. 1982). Triaxial galaxies have principal axes a, 
b, and c of lengths A, B, and C such that A>B>C, 
where A=B for oblate and B=C for prolate ellipticals. 
Many attempts have been made to find the orienta- 
tion of the three principal axes with respect to the 
apparent projection of triaxial galaxies (van Albada 
et al. 1982; Steiman-Cameron, Kormendy & Durisen 
1992 and references therein). The first assumption 
is that the radio axis is aligned with a principal axis 
of the system. It is then commonly accepted that 
the dust lane will lie perpendicular to the radio axis 
for the following reasons: firstly, there is a strong as- 
sociation between galaxies with radio jets and those 
with dust lanes, and hence both features are likely 
to be produced by a related mechanism such as the 
gas layer fuelling the central accretion disk. Secondly, 
surveys of galaxies with dust lanes and radio jets have 
shown a strong correlation between the jet axis and 
rotation axis of the dust lane disk (Kotanyi & Ekers 
1979, Mollenhoff et al. 1992). 

Dust lanes and cylindrical gas disks are thought to 
be the result of mergers, as they are common features 
of galaxies with boxy bulges (eg: NGC7332, Fisher, 
Illingworth & Franx 1994), and are often dynamically 
decoupled from the stellar velocity field (Bertola, Bu- 
son & Zeilinger 1992). Steiman-Cameron & Durisen 
(1982) and Steiman-Cameron et al. (1992) modelled 
gas in-falling into a triaxial potential from several 
merger scenarios. They found that the accreted gas 
from the merger of a gas rich galaxy and an elliptical 
would settle into one of the principal planes of the 



resulting galaxy. The initial orientation of the angu- 
lar momentum vector of the accreted gas, relative to 
the principal planes, will determine into which plane 
it settled through differential precession and dissipa- 
tion. The angular momentum of the accreted gas can 
be parallel, anti-parallel or perpendicular to the stel- 
lar rotation. If the figure rotation axis lies in the 
plane of the dust lane, the result is tumbling about 
this plane. Furthermore, if the gas and dust are from 
separate merger or interaction events, then the cylin- 
drical gas disk rotation can be perpendicular, parallel 
or anti-parallel to the dust lane and radio axis. 

Attempts to model stable orbits in triaxial poten- 
tials have shown that the only stable axes of rotation 
are the long and the short axes, a and c respectively. 
Rotation about the intermediate axis is unstable (van 
Albada et al. 1982). If the axis of rotation of the 
dust disk (and hence the radio axis) is the short axis, 
then a stable polar gas disk is therefore confined to 
be orbiting the long axis, and vice versa (Steiman- 
Cameron & Durisen 1982, Merritt & de Zeeuw 1983). 
The figure rotation may be parallel or perpendicular 
to the radio axis but cannot occupy the intermediate 
axis. If the figure rotation is perpendicular to the ra- 
dio axis, the triaxial figure is tumbling. Furthermore, 
if the surface figure is tumbling, it must do so much 
slower than the orbital period of gas in the polar disk, 
so the disk will not become unstable. The apparent 
stability of polar gas disks implies that these galaxies 
are triaxial and that the surface figure rotates much 
slower than the polar gas disk. 

In order to find stable configurations for rotating 
triaxial systems, it will be assumed that the figure 
rotation is in the same sense as the stellar rotation. 
This applies if the mass distribution follows the light 
distribution. If so, then the four possible stable con- 
figurations are given in Table ^ (Merritt & de Zeeuw 
1983 , van Albada et al. 1982). The warped dust 
lane in ESO248-G10 is aligned closest to the major 
axis (Saripalh et al. 1994). This classifies ES0248- 
GIO as configuration 3 in Table |[ Hence the radio 
axis is the short axis and the figure rotation is about 
the long axis. While the actual gas rotation position 
angle has not been measured, the spectrum showed 
fast rotation along the slit at a p. a. of 40°. The radio 
axis is only 11° away at a p. a. of 51°. Therefore the 
axis about which the gas orbits is likely to be in, or 
close to the plane of the dust lane, perpendicular to 
the radio axis. If the gas is from a recent merger it 
may not yet have settled into a principal plane. If 
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the gas orbit has settled and is stable, then the rota- 
tion axis is in the plane of the dust disk and confined 
to the long axis. In that case, the stars and gas are 
tumbling about the long axis with the figure rotation 
much slower than that of the gas. 

It has been argued that warped dust disks are 
transient and have not yet settled into their rotation 
plane. However, van Albada et al. (1982) showed 
that slow figure rotation or tumbling about an axis 
in the plane of the dust disk can result in a stable 
warped disk due to the Coriolis force. The result- 
ing warp increases with radius. Furthermore, if the 
warped dust lane is along the major axis, then the 
rotation of the warp is prograde to the figure rotation 
when viewed down the figure rotation axis (Merritt 
k De Zeeuw 1983, Kormendy & Djorgovski 1989). 
While the gas and stars are both rotating about the 
long axis, the sense of the stellar rotation is unknown. 
It would be plausible for both to rotate in the same 
sense, although it is not impossible to find counter- 
rotation about the same axis (eg. NGC7332, Fisher 
et al. 1994). If the gas and stars in ESO248-G10 are 
rotating in the same sense, then its figure rotation 
is away from the observer to the north-east. A sta- 
ble system would therefore imply that the south-east 
dust lane warp is rotating away from the observer. 
Furthermore, the optical images reveal the strongest 
extinction from the dust lane is north of the nucleus 
(Figure ||(a)) and hence the dust lane tilts away from 
the sky plane to the north-east. As the warp in the 
dust disk is expected to increase with radius, this con- 
fines the long (figure rotation) axis to have a position 
angle greater than 107° (from the inner dust lane). 
Furthermore, the projected long axis is likely to be 
close to perpendicular (p. a. = 130°) to the slit orien- 
tation as the measured velocity is so fast. Figure ^ 
shows the orientation of the observed axes. The long 
figure rotation axis is shown at a nominal p. a. of 119°, 
midway between the above values. 

4.1.2. A Model for the Orientation of ESO248-G10 

The radio axis is expected to align with t he long 
or the short principal axis (see section [l.l.l] ). How- 
ever, due to orientation effects, the radio axis does 
not appear coincident with the apparent major or 
minor axes. In such a triaxial system, to determine 
the orientation of the principal axes and their rela- 
tive lengths requires five parameters. Firstly, we de- 
fine V' as the apparent angle between the minor axis 
and the projection of the short principal axis on the 



sky. Two more parameters define the principal axes 
lengths such that ( ^ B/A and £, = C/A. The re- 
maining viewing angles are given by the Euler angles 
(6',0). We define two coordinate systems in Figure ^ 
(a,b,c) are the axes associated with a triaxial ellip- 
soid, while (x,y,z) define the sky plane with z along 
our line of sight. 6 is the angle of axis c to the line 
of sight. The angle from axis a to the intersection of 
the c— z and a— b planes defines (j). 

The relationships among these five parameters have 
been extensively discussed in the literature since Stark 
(1977) derived the apparent brightness distribution of 
a set of similar triaxial ellipsoids (see also van Albada 
et al. 1982; de Zeeuw & Franx 1989). We extend 
the notation and equations of Hui et al. (1995) to the 
case of ESO248-G10, where the radio axis is along the 
short axis, and ip is 34° (Saripalli et al. 1994). It can 
be shown that 



^ = v/0.404Lcot6'cosec6' - Mcot26' + iVcosec2£ 



where 



M 
N 
L 



(C^ sin^i/) -I- cos^f; 
(sin^c/) -I- C^cos^f; 
sin20(C^ - 1) 



(1) 

(2) 
(3) 
(4) 



and ip is the apparent angle of the short axis to the 
minor axis as defined above, when 

Mcos^e +(^sm^e - N + 2.475 L cose < 0. (5) 

We now set about finding values of 9, 4>, and 
C that satisfy equations |^— ^. To constrain these 
four parameters, several assumptions will be made. 
As the radio structure appears so large (Saripalli et 
al. 1994), we assume that the radio axis (and hence 
the short axis) is close to the sky plane, say, 9 > 
60°. However, with the radio axis perpendicular to 
the dust lane (see section |]), the tilt of the dust lane 
suggests the radio axis is projected away from the 
sky plane to the north-east, and hence 9 < 90°. We 
therefore adopt 60° < 6* < 80°. There are then 
two constraints on cj). Firstly, for a given 9, 4> must be 
chosen to give the correct apparent angle between the 
figure rotatio n (lon g) and radio (short) axes, called (3. 
From section 4.1.1, the figure rotation axis has p. a. > 
107° and the radio axis has p. a. = 51°, therefore (3 > 
56° . Secondly, as 6 > c, equation |l| requires that for 
any given 0, (f) is limited by C > Figure ^ shows 
the set of curves for which C > ^. For any chosen 
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Table 3: Four possible configurations for rotating triaxial systems. 



Configuration 


Figure Rotation 
(Stellar Rotation) 


Rotation of Dust Disk 
(Radio Axis) 


Form of Dust Lane and 
App. Major/Minor Axis 


Example 


1 


short axis 


short axis 


flat - major 


M31 


2 


short axis 


long axis 


warped - minor 


CenA, M84? 


3 


long axis 


short axis 


warped - major 


M84?, ESO248-G10 


4 


long axis 


long axis 


flat - minor 






ghgjmQS 20'"59" 20'"5B" 20=57= 30"56" 80'°55" 



Right Ascension (J2000) 



Fig. 6. — Orientation of observed axes overlaid on the B— band image and contours. 
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Fig. 7. — Viewing angles 6 and <j) represent the angles between the principal axes, a, b, and c, and the viewing 
axes, X, y, and z, as shown. The z axis is chosen to be along the line of sight, putting the x and y axes in the sky 
plane. ■0 is apparent angle between the projection of the minor axis (y) and the short principal axis, c, on the sky. 



9, there is a maximum possible (j). This restricts the 
range of values to those shaded in Figure ||. Modelling 
has shown that choosing the maximum (j) value puts 
the figure rotation axis approximately perpendicular 
to the slit at p.a. of 132 ± 2°, and /? = 81 ± 2°. We 
can therefore constrain the figure rotation (long) axis 
to be between 107° and 132°; we (arbitrarily) adopt 
the midpoint of 119°, which gives (3 = 68°. 

The dust lane lies in the plane perpendicular to 
the radio axis (see section [4.1.l[ ). 9 can be chosen 
to reproduce the apparent tilt of the dust lane in the 
perpendicular plane. The warp in the dust lane makes 
the tilt difficult to judge, but within the range 60° < 
6* < 80°, modeUing suggests 6* = 65°. For (3 = 68°, cj) 
is found to be 19 ± 1°. 

Once the angles are decided, the dependence of the 
axis ratios is shown in Figure ^. The intrinsic ratio 
C/A > 1/2 is implied by the small number of radio 
galaxies flatter than E5 (Hummel 1980, referenced in 
van Albada et al. 1982). Therefore 0.5 < ^ < 1 and 
C > The choice of ( and ^ affect the apparent 
photometric major/minor axes ratio. Saripalli et al 
(1994) found this axis ratio to be 1.18. Modelling 



an ellipsoid with C = 0.75 ± .05 and C = 0.69 ± .06 
at the chosen angles was found to reproduce this 
major /minor axis ratio, giving 1.18 ± 0.01. 

Figure |l^ shows our adopted model for ES0248- 
GIO with V = 34°, 9 = 65°, cj) = 19°, C = 0.75, 
and ^ = 0.69. While not the only possible model, it 
does represent a self consistent choice of angles and it 
reproduces the observed features. 

4.2. Nature of the Nuclear Source and Cir- 
cumnuclear Structure 

4-2.1. Reddening Processes 

The extinction images in Figure |^(d) reveal sig- 
nificant reddening around the nucleus. Thermal gas 
emission, hot dust emission and dust extinction can 
all affect the galaxy colours. To distinguish between 
these processes, the two-colour pixel plot in Fig- 
ure |ll| was formed from the values of 31 x 31 pixels 
(15.5 X 15.5arcsec) around the nucleus; the typical 
colours of E/SO galaxies and QSOs are shown on the 
figure. Also shown are vectors representing 5 magni- 
tudes of visual extinction, hot dust at temperatures 
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P =56 degrees 



Fig. 8. — ^ vs (j). The shaded region represents the range of values allowed by equation]^ and the constraint C > 
for the chosen range of 9. 



1.0 



5 



major 



major 
minor 




Fig. 9. — £, C for the chosen values of 6, ip, and 
and C as shown. 



The apparent major-to-minor axis ratio is dependent on ^ 
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of 800 and lOOOK mixed with a late-type stellar pop- 
ulation, and thermal gas comprised of hydrogen and 
helium continuum and emission lines mixed with late- 
type stars (Alonso-Herrero et al. 1998). Aperture 
photometry colours for a range of increasing apertures 
are also plotted in Figure 

Notice that in Figure |ll| both the individual points 
coincident with the nucleus and the smallest aperture 
photometry point are bluer than points at larger radii. 
With increasing radius, the colours initially become 
redder along an almost horizontal vector. Beyond a 
radius of 4 arcsec, the points then follow the direction 
of the hot gas vector. Two scenarios could explain 
this pattern. The colours may be the sum of a hot 
dust and a hot gas vector. As the radius increases, 
the contribution of the hot dust and hot gas both 
grow until, at 5 arcsec, the hot gas processes become 
dominant and there is no further contribution from 
hot dust. Similarly, dust extinction, not emission, 
could contribute to reddening out to 5 arcsec. From 
the two-colour pixel plot alone, it is not possible to 
distinguish between these two scenarios and we now 
examine in turn the evidence for dust extinction and 
hot dust emission. 



{a)Dust extinction: This is plausible since the optical 
images clearly show a dust lane. The thickest portion 
of the dust lane does not appear to obscure our view of 
the nucleus, as the tilt results in the foreground dust 
passing above the nucleus. Two observations support 
the interpretation of the two-colour pixel plot as indi- 
cating dust extinction plus thermal gas. Firstly, the 
nuclear pixels do not show as much dust extinction 
as those at larger radius. Secondly, the aperture pho- 
tometry colours follow a dust extinction plus thermal 
gas emission vector until the 4 arcsec aperture point. 
This corresponds to the distance to which the dust 
lane extends from the nucleus in the direction where 
the dust lane is closest to the nucleus in the opti- 
cal images. Beyond 4 arcsec, the aperture photome- 
try points follow the thermal gas vector and show no 
further contribution from dust extinction. Further- 
more, the individual points within the dotted circle 
in Figure ^l], which have the highest reddening, are 
the pixels on the dust lane; dust extinction provides 
a logical explanation for these points. If so, the ex- 
tinction to the dust lane pixels would be around 13 
visual magnitudes. 

{h)Hot dust emission: The second explanation of the 
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0.5 1 1.5 2 2.5 3 3.5 

H-K 



Fig. 11. — J — H YS H — Kn for ESO248-G10. Typical elliptical galaxy and QSO colours are marked by the two 
ellipses. An extinction vector of Ay = 5 is shown along with hot dust vectors for 800K and lOOOK. The colours 
of 31 X 31 pixels (15.5 x 15.5arcs(;c) centred on the nucleus arc shown by the small dots. The crosses are the 30 
pixels closest to the nucleus. Aperture photometry colours marked by large dots represent apertures of 0.5, 1, 1.5, 
2, 2.5, 3, 3.5, 4, 4.5, 5, 6, 7, 8, 9, lOarcsec from left to right. For clarity only some of these aperture radii have 
been labelled. Points within the dotted circle all lie on the dust lane. The dashed line is parallel to the thermal 
gas vector. Dotted lines indicate a themal gas vector which could combine with the dust extinction or hot dust 
vectors to give the observed aperture photometry points. 
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two-colour pixel plot involves hot dust emission as 
the prominent reddening source out to 4 arcsec. Four 
mechanisms may contribute to hot dust emission: 
(i)transient heating of small grains by reprocessed X- 
rays, (ii)scattered nuclear radiation, (iii)heating by 
UV photons from the central source, or (iv)in situ 
heating of dust by hot stars. We now examine each 
mechanism in turn. The first option cannot be ruled 
out in this models of X-ray heating of small 

grains are inconclusive. The processes of grain de- 
struction and heating are not well enough understood 
to produce consistent models of this process. For the 
second mechanism, hot dust emission close to the cen- 
tral source can be scattered into the line of sight by 
clouds further from the nucleus. A polarisation image 
would be required to prove scattering. However, mod- 
els of similar galaxies have shown that essentially all 
the emission would need to be scattered towards the 
line of sight (Alonso-Herrero et al. 1998). Hence, this 
scattering scenario is unlikely. The third mechanism 
for hot dust emission is heating by UV photons from 
the central source which will reradiate in the near in- 
frared. UV heating has been considered to be a major 
source of near-infrared emission from regions close to 
the nuclei of Seyfert galaxies (Alonso-Herrero, Ward 
& Kotilainen 1996; Kotilainen et al 1992; Bryant & 
Hunstead 1999). If the same process is responsible 
for extended thermal dust, then it must be shown 
that UV photons from the central source can heat 
the dust to the observed temperatures out to the ra- 
dius required. For a dust optical depth of Tuy, a total 
UV luminosity of Ijuv (in units of 10^^ erg s~^) would 
heat the dust to a temperature T (K) out to a radius 
r (pc) given by. 



T ■ 



(6) 



(Barvainis 1987). The near- infrared colours indicate 
that if hot dust is the reddening mechanism, then this 
dust extends to ^ 6.5 kpc 4 arcsec) before hot gas 
becomes the dominant mechanism. The typical tem- 
perature for such hot dust is around 800K. A higher 
temperature would produce a steeper vector on the 
2-colour plot. At 1500K the dust would be evapo- 
rated by the UV radiation (Barvainis 1987). There- 
fore, to heat dust to a temperature of 800K out to 
a radius of 6.5 kpc would require a UV luminosity of 
1.5 X lO^^ergs"^, three orders of magnitude greater 
than the most luminous quasars. Therefore, it is very 
unlikely that UV radiation from the central source 
could be responsible for heating dust to such a ra- 



dius. The fourth origin of hot dust emission requires 
groups of hot stars to provide local heating of the dust 
far from the central source. Massive star formation in 
the galaxy halo produces UV emission which can be 
almost entirely converted to thermal dust radiation 
(Alonso-Herrero et al. 1998). Not only can star form- 
ing regions heat the dust, they can also heat the gas. 
As thermal gas emission is an important processes 
further from the nucleus, the evidence for starburst 
regions will be discussed at length in the next section. 
The star forming regions will be shown to extend to 
the radius of the proposed hot dust emission. 

4-2.2. Distribution of AGN and Starburst Processes 

Table | compares three diagnostic emission line ra- 
tios measured in 3 pixels (2.3 arcsec) centred on the 
nucleus with the expected ratios for starbursts and 
AGNs (Radovich & RafaneUi 1996, Hill et al. 1999). 
For each diagnostic, the nucleus is clearly identified 
with AGN emission. 

In Figure |l|(d), horizontal lines represent five re- 
gions where pixels were extracted to make spectra. 
The corresponding spectra and the [N ii]A6584/Hq; 
flux ratios are shown below the rotation curve. For 
radii less than ^ 5 arcsec from the nucleus, the ra- 
tio is > 1.2 as required for AGNs. However, beyond 
this radius, the ratios drop to the range expected 
for starburst/composite regions (< 0.8). Therefore, 
star formation regions are found beyond the radius at 
which the slow rotation 'dip' (see section 3.2.2) starts 
on the rotation curve. It is notable that whereas the 
slow rotation 'dip' is only seen to the north-east of the 
nucleus, the evidence of star formation begins at the 
same radius to both the north-east and south-west. 

4.2.3. North-east 'Dip' in Rotation Curves 

Rotation curve features similar to the north-east 
velocity 'dip' have been observed in other galaxies. 
In each case, however, the explanation for the dip 
does not appear relevant to ESO248-G10. Ionised 
gas rotation curves with a steep gradient followed by 
a dip then a shallower gradient have in some cases 
(eg: NGC253 Arnaboldi et al. 1995) been attributed 
to a nuclear ring. Capaccioli & Longo (1994) suggest 
that stellar rotation curves with a sharp rise up to 
a maximum followed by a dip then a shallower rise, 
are typical of lenticular galaxies (eg: NGC1553, an 
SO galaxy) , in which the dip marks the transition be- 
tween disk-dominated and bulge-dominated parts of 
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Table 4: Measured emission line ratios for three diagnostic emission line pairs, 

measured at the nucleus. 

Line Ratio Measured value Value for pure starburst Value for AGN 

H« — — 

[S II] A 6716+6731 0.73 ± 0.2 < 0.25 > 0.32 

Ha — — 

0.18 ±0.05 < 0.05 > 0.1 

Ha — — 



the galaxy. However, although the ESO248-G10 ro- 
tation curve shows a dip to the north-east, the lack 
of a corresponding dip to the south-west appears to 
rule out a nuclear ring or lenticular galaxy features. 
We now examine alternative models which include a 
second velocity system to the north-east. 

4.3. Possible Models of the Host Galaxy 

A common source of extended emission line gas in 
ellipticals is capture from a merger or interaction with 
a gas rich galaxy (see section [4.1. l[ ). In Figure [T^(a) 
the R image has been divided by the B, to give an 
B — R image in which darker regions correspond to 
redder colours. The nucleus has B — R — 2.5 and is 
marked by a cross. Approximately 10 arcsec to the 
east of the nucleus is a bluer patch with B — R = 1.3, 
first identified by Saripalli et al. (1994). This could 
be a merging companion or a region of blue stars. The 
blue patch appears to be extended towards the west. 
J— band contours overlaid on the B — R image appear 
unaffected by the dust lane to the north-west of the 
nucleus. It is therefore likely that the distortion of 
the J— contours to the east, whilst coincident with the 
dust lane, is due to the extended emission associated 
with the blue patch rather than dust. 

If the blue region in ESO248-G10 can be identi- 
fied with an interacting companion, then its orbit will 
be settling into the principal plane closest to the an- 
gular momen tum vector of the merging galaxy (see 
section 4.1.1 ). In Figure |l^(b), the p. a. of the blue 
patch from the galaxy nucleus is between the figure 
rotation axis and the radio axis. In fact, it could 
be closer to either. If the merging galaxy is in the 
long/intermediate plane, and settling into an orbit 
about the radio axis, then the current merger could 
be responsible for the warped dust lane. While that 



hypothesis can not be ruled out, it would not support 
the distribution of starburst regions or the rotation 
curve features observed. 

An alternative picture involves the merging galaxy 
being located in the short /intermediate plane and 
hence settling into rotation with the figure and gas 
orbits. The dashed lines in Figure |l^(b) are two or- 
bits about the figure rotation axis, projected onto the 
sky plane. Their radii were chosen to lie at either 
end of the putative merging galaxy. In Figure |^(c), a 
region of slower velocity was evident from a radius 
of 5 — 15 arcsec to the north-east, with a local 
minimum at a radius of ~ 8 arcsec from the nucleus. 
A white line along the slit direction in Figure [T^(b) 
marks the radii of the 'dip', with the cross line at the 
local minimum velocity. If gas is being accreted from 
the blue galaxy into an orbit about the figure axis, 
the radii of this gas when rotated to the apparent slit 
angle, corresponds well with the slower gas region. 
However, it is likely that the merging galaxy would 
have an orbit close to the one shown, but not yet 
settled into the short/intermediate plane. Neverthe- 
less, if the blue region is a merging companion, then 
it could be the second velocity system affecting the 
rotation curve. This picture suggests the inner radius 
of the proposed merging galaxy's orbit is currently 
^ 5 arcsec. Hence the radius at which thermal gas is 
accreted into the orbit should be the same either side 
of the nucleus. Newly accreted gas has not yet merged 
with the orbiting gas and could lead to a second veloc- 
ity system near the merging galaxy. Therefore only 
the north-east of the rotation curve is affected by a 
second velocity system. The gas is heated by star- 
forming regions triggered by the merger. The detec- 
tion of starburst emission line ratios should therefore 
be symmetric about the nucleus. This symmetry was 
indeed observed in Figure |l](d) , where starburst line 
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Fig. 12. — {&)B — R image overlaid with J-band contours. The nucleus is marked by the cross. Darker regions 
represent redder colours. There is clearly a blue region to the east of the nucleus. Contour levels are 17.9, 18.1, 
18.4, 18.7, 19.1, 19.3, 19.7, 19.9, 20.1, 20.4, 20.9, 21.2, 21.6 magpixels^^ The J-band has distorted contours along 
the westward extension of the blue patch to the west. This is unlikely to be due to the dust lane as it is not 
apparent in the contours on the west side of the galaxy, (b) B — R image with axes marked. The white line 
along the slit position angle, defines the radii (5-15 arcsec) of the slower velocity 'dip' from the rotation curve. The 
local minimum velocity in the 'dip' is marked by the perpendicular dash across the white line. Two orbits about 
the figure rotation axis are shown by the dashed lines, with radii chosen to lie at either end of the blue feature. 



Rotation in this sense would map the blue feature onto the region of slower velocity (see section |4.3| ) 



ratios became dominant over AGN processes beyond 
^ Sarcsec on both sides of the nucleus. 

An alternative explanation could involve the gas 
originating from a previous merger, and the blue re- 
gion as a cluster of stars. This, however, would leave 
no evidence in any of the images, for the source of the 
second velocity system. Interpretation of the blue re- 
gion as merging galaxy provides a self-consistent ex- 
planation of the observed gas rotation, starburst dis- 
tribution, and optical colours. 

5. FUTURE WORK 

Optical images in better seeing would help to re- 
solve the dust lane features and the proposed merging 
galaxy. Observations to test the interpretation of the 
blue region as a merging galaxy are currently planned, 
including spectral mapping to measure the dynamics 
and distribution of gas across the galaxy. The re- 
sulting spectra across the dust lane may confirm the 
sense of the dust rotation. Furthermore, the rotation 
information across the galaxy would refine the axis 
of gas rotation and our model for the orientation of 
ESO248-G10. 

The radio axis position angle was assumed to be 
that of the radio jets far from the nucleus as mapped 
by Saripalli et al. (1994). However, there is evidence 
that the jet direction may be changing systematically 
with time, possibly indicating precession, and hence 
may not be at exactly the same position angle as fur- 
ther out. Therefore, a detailed radio map of the core 
regions at arcsec resolution could constrain the local 
radio jet position angle and hence further constrain 
the angles found from modelling the system. 

it was shown that there is no [Fe ii] 



In section 



3.2.1 



within 1.6 arcsec of the nucleus, where AGN processes 
dominate. This may indicate that if iron is present 
and not excited by the AGN, then in general, iron ex- 
citation may be more likely from starburst processes: 
in cooling tails of SNRs or photoionizcd by young OB 
stars. An [Fe ii] image and a higher resolution radio 
image of ESO248-G10 could be used to test the ex- 
citation mechanism for [Fell] in active galaxies. The 
emission should be more prominent at a radius > 5 
arcsec (see section 4.2.2 ) if it is due to starburst pro- 
cesses. Furthermore, a correlation between the line 
image and the radio image would be expected for ex- 
citation in cooling tails of SNRs. 



6. CONCLUSION 

ESO248-G10 has been modelled as a triaxial ellip- 
soid with axes ratios of B/A = 0.75 and C/A = 0.69, 
oriented with angles ip = 34°, 9 = 65°, and (p — 19°. 
This is not an exclusive solution but one that is con- 
sistent with the observed apparent axes and rotation 
angles. The core has been identified as an AGN, sur- 
rounded by either hot dust emission or dust extinc- 
tion. Hot gas has been found to play an increasing 
role beyond a radius of 4 arcsec (6.5 kpc). An optical 
spectrum has revealed fast gas rotation out to a ra- 
dius of 20 arcsec (^ 33 kpc). A slower velocity region 
extending from a radius of 5 to 15 arcsec to the north- 
east has been identified as a second velocity feature. 
A model is proposed in which this feature is a merging 
gas-rich galaxy, which is inducing star formation as it 
settles into an orbit about the figure rotation axis. 
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